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Abstract 

We investigate a first-order phase transition from hadronic matter to antikaon condensed matter 
during the cooling stage of protoneutron stars. The phase transition proceeds through the thermal 
nucleation of antikaon condensed matter. In this connection we study the effect of shear viscosity 
on the thermal nucleation rate of droplets of antikaon condensed matter. Here we adopt the same 
equation of state for the calculation of shear viscosity and thermal nucleation time. We compute the 
shear viscosity of neutron star matter composed of neutrons, protons, electrons and muons using 
the relativistic mean field model. The prefactor in the nucleation rate which includes the shear 
viscosity, is enhanced by several orders of magnitude compared with the approximation of earlier 
calculations. Consequently the thermal nucleation time in the approximation overestimates our 
result. Further the thermal nucleation of an antikaon droplet might be possible in our case for 
surface tension smaller than 20 MeV fm~^. 

PACS numbers: 97.60.Jd, 26.60.-c,52.25.Fi,64.60.Q- 
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I. INTRODUCTION 



Antikaon {K~ meson) condensation in dense baryonic matter formed in heavy ion colli- 
sions as well as in neutron stars was first proposed by Kaplan and Nelson There was lots 
of interest in the study of antikaon condensation in neutron stars {2 3] after their work. A 
first-order phase transition from hadronic matter to antikaon condensed matter was inves- 
tigated in several cases using relativistic field theoretical models 0, [ill, isl, Q . The phase 
transition was either studied using Maxwell construction or governed by Gibbs' rules for 
phase equilibrium along with global baryon number and charge conservation l^]. In those 
cases the focus was on the equation of state (EoS) and neutron star structure as well as 
critical temperature of antikaon condensation. 

A first-order phase transition may proceed through the nucleation of droplets of the new 
phase. The formation of droplets of exotic matter such as antikaon condensed matter and 



quark matter, could be possible in neutron stars w^ 
temperature of ~ 10 MeV and is deleptonised [18 



to quantum and thermal nucleation mechanisms 



len the protoneutron star cools down to a 
231 . This nucleation process could be due 
22 1 . The thermal nucleation of antikaon 



condensed matter was already studied using the homogeneous nucleation theory of Langer 



18 



19 



24| . Recently it has been shown in the context of nucleation of quark matter in 



hot and neutrino-free neutron stars that the thermal nucleation is more efficient than the 
quantum nucleation process at higher temperatur es [22 1. 



The homogeneous nucleation theory of Langer 24, l25|] is applicable close to a first-order 



phase transition. The hadronic matter becomes metastable near the phase transition point 
when there is a sudden change in state variables. In this case thermal and quantum fiuctua- 
tions play important roles in the metastable phase. Small ranged and localised fiuctuations 
in state variables of the metastable hadronic matter might lead to the nucleation of droplets 
of stable antikaon condensed matter. Those droplets of the stable phase having radii larger 
than a critical radius will survive and grow. A droplet may grow beyond the critical size if 
the latent heat is transported from the surface of the droplet into the metastable phase. It 
was argued that the heat transportation could be achieved through the thermal dissipation 

The influence of thermal conductivity and shear viscosity on 



and viscous damping 



the thermal nucleation time was studied in a first-order phase transition from hadronic to 
quark matter 



22 



26| . However there is no such calculation for the thermal nucleation of 
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antikaon condensed matter. 

We are motivated to study the effect of shear viscosity on the thermal nucleation rate 
of droplets of antikaon condensed matter in this work. Shear viscosity of pure neutron and 



neutron star matter has been calculated by several groups [28H33|. Recently we have in- 
vestigated the shear viscosity in antikaon condensed matter 3J]. Though we considered a 
first-order antikaon condensation, we did not take into account the nucleation process. We 
performed the calculation of shear viscosity in neutron star matter using the EoS derived 
from relativistic field theoretical models. Now the question is how this shear viscosity of 
nucleonic phase may impact the nucleation of antikaon condensed phase. The earlier calcu- 
ation of the nucleation of quark matter adopted a parametrised form of the shear viscosity 



22| . In this calculation we adopt the same EoS for the computation of thermal nucleation 
time of the antikaon condensed phase and shear viscosity of nuclear matter. 

We organise the paper in the following way. We describe models for homogeneous nucle- 
ation, shear viscosity and EoS in Sec. II. Results of this calculation are discussed in Sec. 
III. Sec. IV gives the summary and conclusions. 

II. FORMALISM 

We consider a first-order phase transition from the charge neutral and beta-equilibrated 
nuclear matter to K~ condensed matter in a hot neutron star after the emission of trapped 
neutrinos. Droplets of antikaon condensed phase are formed in the metastable nuclear matter 
due to thermal fluctuations. Droplets of antikaon condensed matter above a critical size are 
of interest because those critical droplets will drive the phase transition. We adopt the 



homogeneous nucleation formalism of Langer to calculate the thermal nucleation rate 



In this formalism, the thermal nucleation per unit time per unit volume is given by 2^, |25| 



r = V,exv ( ) , (1) 



T 

where AF is the change in the free energy to produce a critical droplet in the metastable 
nuclear matter. The change in free energy of the system due to the formation of a droplet 



is given by 
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23| 



Att 

AF{R) = - — (P^ - P^)P=^ + inaR^ , (2) 
3 
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where R is the radius of the droplet, a is surface tension of the interface separating two phases 
and and are the pressure in nuclear and antikaon condensed phases respectively. We 
obtain the critical radius of the droplet from the maximum of AF{R) i.e. SrAF = and it 
is 

Rc = T,J^ . (3) 



The EoS enters into this calculation through the difference in pressures of two phases. 



The pref actor in Eq. (1) is factorised in the following way 



25 



27| 



2tt 



-fin 



(4) 



The statistical prefactor (Qq) gives the available phase space around the saddle point at Rc 
during the passage of the droplet through it. The expression for the statistical prefactor is 

Here ^ is the correlation length for kaons. This correlation length is considered to be the 
width of the interface between nuclear and antikaon condensed matter. Next we focus on 
the dynamical prefactor k. This prefactor is responsible for the initial exponential growth 

3,H 



rate of a critical droplet. It is given by 

2(7 



Rl,{Awy 



AT + 2(-r/ + C) 



(6) 



Here Aw = wk — wat is the enthalpy difference of two phases whereas A is thermal conduc- 
tivity and 7] and C, are the shear and bulk viscosities of nuclear matter. It was shown that 



the latent heat would 
and viscous damping 



DC taken away from the surface of the droplet due thermal dissipation 
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Finally the thermal nucleation time {Tth) in the interior of neutron stars is calculated as 



(7) 



where the volume V = Att/SR'^^^. We assume that pressure and temperature are constant 
within this volume in the core. It is evident from Eq. (7) that the thermal nucleation 
time is inversely proportional to transport coefficients. In connection with the quark matter 
nucleation, it was noted that the contribution of thermal conductivity is not appreciable. 
Further, we know that shear viscosity increases [sj] and bulk viscosity decreases 35 1 as 
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temperature decreases. After the emission of neutrinos a hot neutron star cools down, and 
the thermal nucleation drives the phase transition when the temperature is ~ 10 MeV. In 
this temperature regime, shear viscosity might be an important factor. 

Now we focus on the calculation of shear viscosity in neutron star matter. It was noted 



that the main contributions to the total shear viscosity in neutron star matter came 
electrons, the lightest charged particles, and neutrons, the most abundant particles 



28 



rom 



32 



30, 



331] . Recently the proton shear viscosity, though small, was considered in our calculation 



34| . We calculate shear viscosities for different particle species using coupled Boltzmann 



transport equations 
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32 



34( 1 . The effects of the exchan ge o f transverse plasmons in the 



32 



36| . The knowledge of nucleon 



collisions of charged particles are included in this case 
nucleon scattering cross sections derived in the Dirac-Brueckner appro ach 
for the calculation of neutron and proton shear viscosities 



32 



34 



37 



38| is used 



39j. Contributions of 



neutrons, protons, electrons and muons to the total shear viscosity 



Vtotal = Vn + Vp + Ve + 



are given by 32|, |34 1 



5m* 



(8) 



(9) 



Here Xj is the relaxation time of i-th species as calculated in Ref. [32|, |3J]. Effective mass 
and Fermi momentum of i-th particle species are denoted by m* and pp., respectively. For 
electrons and muons, effective masses are taken as their chemical potentials due to relativistic 
effects. We adopt a relativistic field theoretical model which is described in the following 
paragraphs, for the calculation effective masses and Fermi momenta of neutrons and protons. 

Besides the computation of shear viscosities, the knowledge of the EoS is essential for 
the thermal nucleation rate. We note that the EoS enters into the calculation of thermal 
nucleation rate as given by Eqs. (l)-(3). We adopt relativistic field theoretical models to 
describe the /3 equilibrated matter in nuclear and antikaon condensed phases. Those two 
phases are composed of neutrons, protons, electrons, muons and of K~ mesons only in the 
antikaon condensed phase. Both phases are governed by baryon number conservation and 



charge neutrality conditions 17|. The baryon-baryon interaction mediated by the exchange 
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of 0", uj and p mesons is described by the Lagrangian density 



42| 



vCat = ^ {ilt^d^ -rnB + g^rBcr - gu^Blf,^^ - gpBli^B ■ P^) ipB 

B=n,p 



+ -{d,ad^cr-mla^)-U{a) 



4"/^— ■ 2 
The scalar self-interaction term 



pv 



BQ, 43| 



is 



The effective nucleon mass is given by = — gaB<^, where is the vacuum 
mass. The Lagrangian density for (anti)kaons in the minimal coupling is given by 0, 



(10) 



(11) 



baryon 



Id, 



12| 



K 



D*^KD''K - m*^KK , 



(12) 



where the covariant derivative is = + WuiK^p + igpK^K ■ Pp and the effective mass of 
(anti)kaons is m*^ = itlk — gaKO'. The in-medium energy of K~ mesons is given by 



gujK^O + -gpKP03 



(13) 



We consider the s-wave (p = 0) condensation. The condensation sets in when the chemical 
potential of K~ mesons (/ix- = ^jJr-) is equal to the electron chemical potential i.e. pe = 
Pk-- The critical droplet of antikaon condensed matter is in total phase equilibrium with 
the metastable nuclear matter. The mixed phase is governed by Gibbs' phase rule along with 
global baryon number conservation and charge neutrality 17|. Two phases are in chemical 
equilibrium and the mechanical equilibrium is constrained by Eq. (3). 

We do not consider the variation of the meson fields in the droplet with position. This 
is known as the bulk approximation 19|. We solve equations of motion self-consistently in 
the mean field approximation 40[ and find effective masses and Fermi momenta of baryons. 
Further we obtain the pressure in nuclear matter (P^) and in antikaon condensed matter 
(P^) as given by Ref.(l2|. 



III. RESULTS AND DISCUSSIONS 



We adopt the Glen denning and Moszkowski parameter set known as GMl 4^ for nucleon- 
meson coupling constants which are obtained by reproducing the saturation properties of 
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nuclear matter such as binding energy E/B = —16.3 MeV, baryon density hq = 0.153 
fm~^, asymmetry energy coefficient Oasy = 32.5 MeV, incompressibility K = 300 MeV and 
effective nucleon mass m*^/mN = 0.70. Next we determine kaon-meson vector coupling 
constants using the quark model and isospin counting rule {3, [s^- And the scalar kaon- 
meson coupling constant is obtained from the real part of K~ optical potential depth at 
normal nuclear matter density. 

Heavy ion collision experiments and K~ atomic data indicated an attractive potential 



for antikaons and a repulsive potential for kaons 45N50|. However the strength of antikaon 



optical potential depth is a debatable issue. It was found from the analysis of K atomic 



data that the real part of the antikaon optical 



MeV at normal nuclear matter density 



45 



potential could be as large as Uj^ = —180 ±20 



46| . But theoretical models including chirally 



motivated coupled channel models as well as double pole structure of A (1405) could not 



find such a strongly attractive antikaon optical potential depth |5ll-l54j. In this calculation 
we take an antikaon optical potential depth of = —160 MeV at normal nuclear matter 
density and obtain the kaon-scalar meson coupling constant g^rK = 2.9937. 

Here we study the thermal nucleation of antikaon condensed phase after the emission of 
trapped neutrinos followed by the evolution of a hot neutron star to the cold star. The stellar 
matter is heated during the diffusion of trapped neutrinos leading to the deleptonisation and 
maximum entropy per baryon s = 2. It would take about a few 10^ s for the maximally 
heated star to cool down to a temperature ~ 1 MeV. The threshold for the appearance 
of antikaon condensate would be reduced after the deleptonisation. We demonstrate this 
evolution of a neutrino-free neutron star with a few snap-shots corresponding to entropy 
per baryon s = 0, 1, and 2. We adopt the finite temperature calculation of Ref. jj] in this 
case. We show the temperature as a function of baryon density in nuclear matter for fixed 
s = 1 and 2 in Fig. 1. The temperature increases with baryon density. This shows that 
the temperature varies from a higher value in the core to a smaller value at the surface of a 
neutron star. The equations of state of neutrino free nuclear matter at fixed s = 0, 1, and 
2 are exhibited in Fig. 2. It is evident from Fig. 2 that the temperature of a few tens of 
MeV does not modify the EoS appreciably upto energy densities ~ 1000 MeV fm^'^ relevant 
at the center of a neutron star, with respect to the zero temperature EoS corresponding to 
the s = 0, case because the temperature is much less than the baryon chemical potential. 
Therefore, we consider the zero temperature EoS for the rest of our calculation. 
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The total shear viscosity in nuclear matter is shown as a function of normalised baryon 
density for different temperatures in Fig. 3. The shear viscosity decreases as temperature 
increases. Further, the shear viscosity increases with increasing baryon density. It was 
earlier noted that the temperature dependence of shear viscosities corresponding to different 
particle species was a complicated one and manifested through relaxation times 3J|. This 
temperature de pen dence deviates from the characteristic temperature dependence, of 
a Fermi liquid |32|. Figure 3 shows the total shear viscosity for three temperatures T = 1, 
10, and 30 MeV. We are interested in the temperature range from a few MeV to a few tens 
of MeV in this calculation. 

Next our focus is on the prefactor (Fq) of Eq. (4) which not only involves transport 
properties of nuclear matter such as thermal conductivity and shear and bulk viscosities but 
also depends on the correlation length of kaons and surface tension. The radius of the droplet 
is to be greater than the correlation length for kaons ^. Otherwise many approximations in 



ation 



22 



ength is also taken as the thickness 



26| estimated to be ~ 5 fm 19|]. Here 



this calculation would collapse [19|, l26[ . The corre 
of the interface between nuclear and kaon phases 
we consider antikaon droplets with radii greater than 5 fm. The other important parameter 
in the prefactor is the surface tension. The surface tension between nuclear and kaon phases 
was already calculated by Christiansen and collaborators [55]. They found that the surface 
tension varied from a value of 30 MeV fm~^ at the start of the mixed phase to a value of 
10 MeV fm~^ at the end of the mixed phase. It is evident from the fit of Ref. 55|] that the 
surface tension is sensitive to the EoS. It is to be noted that we have a different parameter 
set and stiffer EoS than those of Ref. |55|. We perform this calculation for a set of values 
of surface tension a = 10, 15 and 20 MeV fm~^. Further, we approximate the difference in 
enthalpy densities (Aw) in Eq. (6) as the difference of energy densities because the pressure 
difference is negligible compared with the difference in energy densities 26 1. 

The prefactor (Fq) is shown as a function of temperature in Fig. 4. The prefactor is 
calculated at a baryon density rif, = 2.3no and surface tension a = 10 MeV fm~^. The 
density chosen here is just above the critical density for the transition. We calculate the 
prefactor for two values of correlation length ^ =2 and 5 fm. Solid lines correspond to 
the calculation of prefactor including the shear viscosity term and neglecting the thermal 
conductivity and bulk viscosity terms in Eq. (6). On the other hand, dashed lines imply 
the prefactor including only the thermal conductivity term. The thermal conductivity is 



8 



calculated according to Ref. 39|]. The prefactor involving the shear viscosity term is a few 
orders of magnitude higher than that of the thermal conductivity for both values of kaon 
correlation length. It is worth mentioning here that the contributions of thermal conductivity 
and bulk viscosity in the prefactor were neglected in an earlier calculation of quark matter 
nucleation 22|]. We observe that the prefactor for ^ =5 fm is substantially reduced from 



that of ^ =2 fm in both cases. In the following paragraphs we discuss the results of our 
calculation performed with the prefactor involving only the shear viscosity term. 



The prefactor of Eg. 



in some calculations 



23 



4) was approximated by according to the dimensional analysis 



261]. This approximation might differ considerably from the actual 



calculation of the prefactor. Fig. 5 shows the prefactor of Eq. (4) including only the 
contribution of shear viscosity along with the prefactor approximated by as a function 
of temperature. It is evident from Figure 5 that the approximated prefactor is several 
orders of magnitude smaller than the prefactor calculated according to Eq. (4) at lower 
temperatures. However, this difference between two calculations is reduced appreciably at 
higher temperatures. It remains to be seen how two calculations of the prefactor influence 
the thermal nucleation time. 

Now we discuss the nucleation time of a critical droplet of antikaon condensed matter 
calculated with the prefactor of Eq. (4) involving only the shear viscosity contribution. The 
critical radius of the droplet is obtained from the relation in Eq. (3). The droplet which 
overcomes the maximum height of the energy barrier of Eq. (1) occurring at the critical 
radius, grows exponentially. Such a droplet triggers the onset of antikaon condensation in 
the interior of a neutron star. For surface tension a = 10 MeV fm~^, we obtain an antikaon 
droplet with critical radius 6.32 fm at baryon density 2.3no. This density is just above 
the critical density for transition from nuclear matter to the droplet of antikaon condensed 
matter where the volume fraction of the antikaon condensed phase is approaching zero. We 
perform the calculation of thermal nucleation of the critical droplet within a volume with 
Rnuc = 100 m in the core of a neutron star as mentioned in Eq. (7). It is assumed that 
density, pressure and temperature are constant within this volume. The thermal nucleation 
time is plotted with temperature in Fig. 6. Results shown here are calculated with the 
kaon correlation length ^ = 5 fm. For the above-mentioned case, the nucleation time of 
the critical droplet drops as temperature increases. For the case of cr = 10 MeV fm~^, 
the nucleation time is 1 s at a temperature of T ~10 MeV. We also estimate thermal 
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nucleation times of critical antikaon droplets for two other values of surface tension a = 15 
and 20 MeV fm~^. It is noted that the size of the critical droplet increases with increasing 
surface tension. Radii of the critical droplets are 9.48 and 12.6 fm corresponding to a = 
15 and 20 MeV fm~^, respectively, at a density 2.3^0. We find a similar trend of the 
nucleation time with temperature in these cases. However, the temperature corresponding 
to a particular nucleation time for example, 1 s, increases as the surface tension increases. 
It is worth mentioning here that the condensate might melt if the temperature is higher 
than the critical temperature IgI. We compare thermal nucleation times corresponding to 
different values of the surface tension with the cooling time of a protoneutron star tcooi ~ a 
few 10^ s. Thermal nucleation of the antikaon condensed phase may be possible when the 
thermal nucleation time is less than the cooling time. It is evident from Fig. 6 that the 
thermal nucleation time is strongly dependent on the surface tension. Thermal nucleation of 
an antikaon droplet is a possibility for surface tension a < 20 MeV fm~^ because in this case 
the thermal nucleation time is less than the cooling time at^ a temperature which might be 
below the critical temperature of antikaon condensation [16]. We repeat the calculation of 
the thermal nucleation time with the correlation length ^ =2 fm and find that the thermal 
nucleation time is slightly higher than that of ^ =5 fm over the whole temperature regime 
considered here. 

We compare the results of thermal nucleation time taking into account the effect of shear 
viscosity in the prefactor with that of the prefactor approximated by T^. These results are 
shown in Fig. 7 for surface tension a = 10 MeV fm~^ and at a density rib = 2.30no . The 
results of the approximation are a few orders of magnitude higher than those of our 
calculation. For example, the thermal nucleation time is 1 s at T ~ 10 MeV in our case, 
whereas it is ~ 10^ s at the same temperature for the approximated case. These results 
demonstrate the importance of including the shear viscosity in the prefactor of Eq. (4) in 
the calculation of thermal nucleation time. 



IV. SUMMARY AND CONCLUSIONS 



We have investigated the first-order phase transition from nuclear matter to antikaon 
condensed matter through the thermal nucleation of a critical droplet of anikaon condensed 
matter. Our main focus in this calculation is the role of shear viscosity in the prefactor and 
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its consequences on the thermal nucleation rate. We adopt the same EoS derived in the 
relativistic mean field model for the calculation of shear viscosity and thermal nucleation 
rate. In this connection, we have constructed critical droplets of antikaon condensed matter 
above the critical density for different values of surface tension. Droplet radii increase with 
increasing surface tension. We obtain thermal nucleation time as a function of temperature 
for a set of values of surface tension and find that the thermal nucleation time is strongly 
dependent on the surface tension. Our results show that thermal nucleation of a critical an- 
tikaon droplet could be possible during the cooling stage of a hot neutron star after neutrinos 
are emitted for a lower value of surface tension a < 20 MeV fm~^. Further, a comparison 
of our results with that of the approximation shows that the approximation overesti- 
mates our results of thermal nucleation time. This comparison highlights the importance of 
shear viscosity in our calculation. 

We have considered the thermal nucleation in a nonrotating neutron star and neglected 
the bulk viscosity in the prefactor. However neutron stars rotate very fast at birth. They 
emit gravitational waves and become r-mode unstable jsSj. The bulk viscosity plays an 
important role in damping the r-mode instability. The bulk viscosity might dominate over 
the shear viscosity above 10 MeV [35]. It would be interesting to investigate the effect of 
bulk viscosity on the thermal nucleation time in a rotating neutron star along with that of 
shear viscosity. 
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FIG. 2. Pressure is plotted with energy density for neutrino-free nuclear matter and fixed entropy 
per baryon s = 0, 1, and 2. 
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FIG. 3. Total shear viscosity in the nuclear matter phase is plotted with normalised baryon density 
for different temperatures. 
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FIG. 4. Pref actor is shown as a function of temperature for different correlation lengths and at a 
fixed density and surface tension. 
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FIG. 5. Same as Fig. 4, but the prefactor including the contribution of shear viscosity is compared 
with that of approximation. 
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FIG. 6. Thermal nucleation time is displayed with temperature for different values of surface 
tension. 
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